VORTEX COLLAPSE FOR THE ^-CRITICAL NONLINEAR SCHRODINGER 

EQUATION 



G. SIMPSON & I. ZWIERS 



Abstract. The focusing cubic nonlinear Schrodingcr equation in two dimensions admits vortex solitons, 
standing wave solutions with spatial structure, Q( m '(r, 9) = e lme fl( m ) (r). In the case of spin m = 1, we 
prove there exists a class of data that collapse with the vortex soliton profile at the log-log rate. This extends 
the work of Merle and Raphael, (the case m = 0,) and suggests that the L 2 mass that may be concentrated 
at a point during generic collapse may be unbounded. Difficulties with m > 2 or when breaking the spin 
symmetry are discussed. 
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1. Introduction 

We consider the L 2 -critical nonlinear Schrodinger equation in two dimensions, 

f iu t + Au + u\u\ 2 = 
(1.1) ^ 

[ u(0,x) = u G H^R 2 ). 

Equation ( |l.l| is locally wellposed for data uq G H 1 , [7 11 . That is, there exists a solution u G 
C ([0, T max ), i? 1 ) and some fixed negative power so lino || m . Therefore, we have 

the classic blowup alternative, 

T max = +oo or, lim ||u(i)|| H i = +oo. 

t-f-Tmai 
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Evolution of uq by equation (1.1 1 preserves the following quantities. 

(1.2) M[u ] = M[u(t)) = \u(t,x)\ 2 dx, (mass) 



(1.3) E[u ] = E[u(t)]= j\V x u(t,x)\ 2 dx-^J\u(t,x)\ 4 dx, (energy) 

(1.4) P[uq] = P[u(t)] = Im I / u(t,x)Vu(t 7 x) dx I . (momentum) 



The associated symmetries of the equation are phase, time translation, and spatial translation. There is a 
Galilean symmetry, 

up Q (t,x) = u(t,x - /3ot)e 1,S z"( x ~^' t \ for any fixed /3 € K 2 , 
and a scaling symmetry, 

u\ Q (t, x) = Ao«(Ao£, Xox) for any fixed Ao > 0. 

The effect of scaling on Sobolev norms is, ||wA ||jj s = AQ S ||ii|| ,j s , for any reasonable s. Note that only the 
critical norm is left invariant. By choosing Ao = ||ii(t)||#i at a fixed time, and using the minimum local 
wellposedness time for unit data in H\ we have the scaling lower bound for the blowup speed, 

u(t) e C (\p,T mayi ),H l ) , with T max maximal, then \\u(t)\\ H i > 1 

V ^ max £ 

Alternatively, the scaling lower bound can be established through energy conservation, [2]. 
Peculiar to the L 2 -critical case, there is also the pseudo-conformal (or lens) symmetry, 

1 / 1 x 

(1.5) v(t.x) = — u — — rrr, — 

\ ) T — t \(T—t) 2 T-t 

which acts on the virial space, {/ G if 1 } PI {|x| 2 / £ L 2 }. In particular, the pseudo-conformal symmetry 
transforms standing wave solutions into blowup solutions with H 1 norm growth jr^. 



1.1. Blowup with Soliton Profile. To find standing wave solutions of equation (1.1), introduce the usual 
ansatz, u(t,x) = e lt Q(x), to derive the profile equation, 

(1.6) AQ-Q + Q\Q\ 2 = 0. 



There is a unique real- valued positive radial solution Q to equation (1.6), as proved by McLeod and Serrin 
16[ This solution we call the soliton, or the ground-state since E(Q) — 0. In this paper we will focus 



on other solutions of equation (1.6), as we discuss in the next section. Weinstein 31 identified the soliton 



as the unique minimizcr of J[f] = — rjri — — among H 1 functions, thereby showing the optimal constant of 



the Gagliardo-Nirenberg inequality, 



\vf L 4 < ||g,| 2 IMfeiHla- 



L' 2 

Note that if M[«o] < ^[Q], the Gagliardo-Nirenberg inequality gives apriori control of the H 1 norm from 
the conservation of energy. That is, there is g lobal wellposedness for data with M[u ] < M[Q]. 

The pseudo-conformal transformation (1.5) applied to the standing wave solution e u Q(x) gives an explicit 
blowup solution with M[uq] — M[Q}. We denote this explicit solution S(t); Merle [17) showed that, up to 
symmetries, it is the only blowup solution with the mass of Q. Bourgain and Wang 111 proved that S(t) is 
stable with respect to perturbations that are exceptionally flat near the central profile. 

More generally, negative energy data in the virial space leads to blowup, as shown by Glassey's virial 
identity [8], 

^ J M 2 K*)I 2 = A j t Im J x • = 16E i u o] 

Ogawa and Tsutsumi 25 later extended the argument to negative energy radial data. 



following earlier work by Coffman [ij in 3D. Kwong 12 extended the result to all H 1 -subcritical nonlinearities. 
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Let us consider B a = {u e H 1 : M[Q] < M[uq] < M[Q] + a}, where a > is some small constant. 
Merle and Raphael [19] proved that there is no solution in B a that blows up as predicted by Glassey's virial 
identitjJ^J They also showed 18 22 that there is an open subset O C B a , including all the negative energy 



data, that lead to blowup in finite time with the log-log rate 

\\u(t)\\ Hl r 



log|log(T-i)| 



T-t 



Raphael 27 proved that all solutions in B a that lead to blowup either belong to O, or blowup with at least 
the H 1 growth rate of S(t). Finally, Merle and Raphael 20 showed that all solutions in B a that blowup 
concentrate exactly the profile Q at a point, in the sense that there are parameters X(t) > 0, j(t) 6 K and 
x(t) £ M 2 such that, 



u{t,x)- W) Q 



x(t) 



A(t) 



-i 7 (t) 



u*(x), 



where the convergence is in L 2 as t — > T" max . Moreover, the residual profile u* identifies the blowup regime, 
with u* H 1 if and only if the solution belonged to O and followed the log-log rate. 



1.2. Vortex Solitons. Vortex solitons are solutions to equation (TL6l of the form Q^(r, 9) = e m9 R^ m \r), 
where R\ m ^ is real-valued and positive. That is, we seek a function Ry 71 ' that satisfies, 



(1.7) 



77! 



R (m) 



0, 



a r i? (m) | r=0 = 0, R (m) (\x\) > 0, R {m) £ H^M. 2 ) n Ux\- l f(x) e L 2 } 



For all m £ Z, Iaia and Warchall M showed there exists a solution to (1.7| and, analogous to the result of 
Kwong [12] in the case m = 0, Mizumachi [23] has shown it is unique. Fibich and Gavish |4j Lemma 12] have 



remarked that the resulting profile Q m is the unique minimizer of </[/] 



\vf\i 2 \f\i 



among H 1 functions 



with spin m. We denote this space by H( m y Some vortex solutions are pictured in Figure jTJ and their radial 
profiles appear in Figure [2] 

This variational characterization gives an optimal Gagliardo-Nirenberg inequality for functions in Hly 

As a consequence, for data u £ H} m \ and L 2 norm less than ||Q( m ) ||^2 there is global wellposedness. As 

a second consequence, Fibich and Gavish [4l Corollary 16] remark that \\Q 



sequence in m. Indeed, Pego and Warchall [26] showed the asymptotic form, 



(m)||2 

'II 2 



is a strictly increasing 



R {m} {r) 



V^sech 



y ^max) 



where r max « ^/2m for m^> 0. Therefore, HQ^H 2 ^ 
to be a good approximatiorj^ve n for small m. 
The linearization of equation near Q( m ' ) is, 



4V3m for large m, which Fibich and Gavish found 



(1.8) 



d t v 



iL (m) [vl where, L^ m \v] = f-A + 1 - \Q 



)(m)|2 



2 Q 



(m) 



Written as a harmonic series, v — X^ez e^ m+ ^' e ,(r), 



(1.9) 
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-A + 1 - |Q (m) | 2 ) e^+Wfj - 2\Q 



)(™)|2 J(m-j)e 



so that it is clear the linear system excites harmonics in pairs. In the case involving only j 

~im8 



< ""*" («i + iv 2 ), we may write — [v] in matrix form as, 



0, that is, 



(1.10) 



o L { r } ~ 






-L { f 






a for which limt_ 







•^Error less than 3% for n 



where, 



:/N 2 |«W| 2 

2, less than 0.4% for m > 5. 



L { l n) = -A - 
L {m} = -A - 



l-3|Q( m )| 2 , 
1 - |Q (m) | 2 . 



0, in constrast to the explicit solution S(t). 




Figure 2. The radial component of some vortex solutions. 



Comparing equations 1 1.9 1 and ( 1.10 1 we see that takes on the form of ( 1.10 1 on all of H 1 in the case of 
spin m = 0. In this important case, Weinstein |32| showed that the generalized nullspace of L has dimension 
8 and is generated by the symmetries. 

In the cases m = 1 and m = 2, the generalized nullspace of L^™) is generated in the same way. However, 
in these cases Pego and Warchall 26 found unstable eigenvalues and additional eigenvalues in the spectral 



gap (all for modes with \ j\ ^ 0). That is, there exists a function p with spin m — 1 such that, 

L { r } (P) = -\y\ 2 Q m , L {m) (\y\ 2 Q m ) = -4AQ m , 

L { l n) (AQ m ) = -2Q m , L {m) (Q m ) = 0, 

where A = 1 + y ■ V denotes the scaling operator. The remaining Jordan chains, generated by X7Q( m \ 
consist of functions with = 1. Instability of vortex profiles is not restricted to the cubic nonlinearity. 



Mizumachi 24 has shown that there are unstable vortex profiles for any power-type nonlinearity strictly 



stronger than linear. 

1.3. Blowup with Vortex Profiles. Any vortex soliton becomes a blowup solution through the pseudo- 
conformal transformation. Study of the asymptotic profile during vortex blowup was initiated by Fibich & 
Gavish [4j, including the variational structure referenced above. Their work includes numerical simulations 
where they found data with mass slightly larger than Q( m ) that blowup at exactly the scaling lower bound 
and with profiles different from the vortex soliton. 

Our main result is that there is a class of solutions with spin m = 1 that blowup with exactly the vortex 
soliton profile and log- log behaviour similar to that established in the case m = 0. 

Theorem 1.1 (Log-log Blowup with Vortex Profile). Assume the Spectral Propertwjs true for spin m. Then 
there exists a class of data T>^ m \ open as a subset of H}y such that for u G V^ m ) the evolution u(t) by 



(1.1) blows up at finite time T max with the Q' m ' profile and log-log rate. That is, fort G [0, T max ) there exist 
continuously variable parameters X(t) > and j(t) G M with the following properties: 
Log-log Blowup Rate:: 

(1.11) lim |Kt)|Un/ - TmaX-t r=C 

K ' t^T m J log|logT max -i| 

Description of the Singularity:: 

(1.12) t Jm^u(t,a;) - e"^ = u*(x) e i 2 (M 2 ). 

We will now discuss the consistency of the self-similar regime discovered by Fibich & Gavish and Theorem 
|1.1| Consider, 

B a<m - {u Q G H\ m) : M[Q(™)] < M[u ] < M[Q (m) ] + a} . 
Then, due to the variational characterization of Q^: 

Theorem 1.2 ("Orbital Stability"). For a > sufficiently small, let v G -B Q , m with, E[v] < a\\v\\ 2 H1 . Then 
there exists A„ > 0, 7^ € R such that, 

\\\ v v{\ v y)e^ -Q {m \y)\\ H , <5{a), 

where 5(a) — > as a — > 0. 

The proof is by means of concentration compactness and the Gagliardo Nirenberg inequality in H} m \ , and 
is not constructive. See 28 Theorem 6] for a clear exposition. The class of data from Theorem |l.l| 



belongs to B a ^ m , and we note that the orbital stability of Theorem 1.2 applies to all data vq G B a ^ m that 
blowup in finite time. Indeed, we conjectur^hat finite time blowup solutions from the class B a ^ m either 
obey the log-log blowup rate or the lower bound, ||«(t)||.ffi > 

(T max t) 



In particular, they present results using uq = 1.02Q' m ' (r, 9) and spin in = 2. The profiles identified, denoted G n 
truncated solutions of equat ion l|2.18| l with Cauchy boundary conditions and an implied value of b, in this case b PS 0.1092. Our 
own truncated solutions of l |2.18| l are very similar. Fibich & Gavish have conveyed by personal communication corresponding 
discoveries for spin m = 1 and data as small as l.OOOOlQ'" 1 '. 

^See Proposition 1 1.1 1 below. 

^We expect the analysis of [27] to apply, and that the proof of Theorem 1 1 . 1 1 may be reformulated to apply to all vq £ B a m 
with \\v{t)\\ H i e L 1 (t e [0,T m ax)), as in [22]. 



Collapse at the square-root rate has also been observed numerically in the case with no spin, [5]. T hese 



examples are an important area of continuing study. It is possible that the threshold a of Theorem 1.2 (and 
hence the applicability of Theorem |1.1[ ) is exceedingly small. 

1.4. Spectral Propety. In order to demonstrate the dynamic claimed in Theorem |1.1[ we will attempt to 
parameterize the solution in terms of the symmetries and a suitable deformation of the profile Q^ m \ In order 
for the finite-dimensional system of parameters to capture the essential dynamics of the solution we require 
two things. First, that the parameter dynamics can be reliably predicted from a finite system of differential 
inequalities. Second, that after removing the central profile from the solution the error e can be estimated 
in terms of those parameter dynamics. 

That the parameter dynamic are stable is an essential feature of the log-log regime. Indeed, Raphael 
showed 27 that the relationship between a particular ratio of parameterfFknd a fixed constant evolves 



according to a Riccati equation, with the log-log dynamic corresponding to the stable branch. To control 
the error e in terms of the dynamics, we will consider the following operator, derived from the linearized 
energy, 



(1.13) H (m \e,e) = (c { ; n) e 1 ,e 1 ) + (c {m 



2 £ 2,£2 



where 



(1.14) 4 m) = -A + 3Q("VvQ (m) , 4 m) = -A + Q("VVQ (m) , 

(1.15) e_ = e lme Re (e~ mie e) , e 2 = e lm0 Im (e" mie e) . 

This decomposition, e = t\ + ie_, is powerful, as it reduces the algebraic structure of the problem in 



to that of the radially symmetric problem in H . For further discussion, see (2.32 1, below. 
We will prove the following for m = 1, 

Proposition 1.1 (Spectral Property). Let e 6 H} m ) Then there exists a universal constant S m such that 

U^^e) >6 m (J \V v e\ 2 + J |e| 2 e-l^ 

(1.16) ^(e u Q^)\(e u kQ (m) ) 2 

+ (e 2 ,AQ (m) ) 2 + (e 2 ,A 2 Q( m )) 2 

In the case of L 2 -critical nonlinearity, no spin, and dimension N = 1, Merle and Raphael [2l[ Appendix 
A] gave an explicit proof of the Spectral Property. In the case of L 2 -critical nonlinearity, no spin, and 
dimensions N = 2,3,4,5, including equation (1.1 1 in the case m = 0, Fibich, Merle and Raphael [6] have 



given a numerical proof that inspires our own proof of Proposition 1 1 . 1 1 in Section [3j Details of our numerical 
methods are provided in Appendix[j5j Code to reproduce our computations is available at http : //www . math 



toronto.edu/simpson/files/vortex_dist.tgz. As stated, the spectral property is false for m = 2,3. 



2. Proof of Log-log Blowup 

In this section, we prove Theorem |1.1| assuming Proposition |1.1| Before decomposing the solution, we 
introduce almost self-similar deformations of the vortex profiles that simulate the effect of symmetries that 

do not belong to H 1 . The standard self-similar ansatz is, u(t, x) — , 1 I , x | e %u ^\ which 

6 ' V ' ' y/2b(T-t)^ b \y/2b(T-t) J 

gives the following equation for the spatial profile, 

(2.17) AQi m) - Q { b m) + l bKQt n) + Q { b m) IQ^ ? = 0- 



Namely the sign of /— = x — do V Eo . Parameter b will be introduced in Section 
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2.1 



We seek solutions with spin to. Remove a quadratic phase, e tm9 p^ m ' ( r ) = QY n> e * , and assume the radial 
profile (r) is real valued. We seek solutions to, 

b 2 



(2.18) 



Ap M _ [ 1 



- -r z P 



= 0, 



lim r- m P^ m \r) ^ 0, lim d r 
r _>0+ r->0+ 

As pointed out by Fibich and Gavish [ij Lemma 8], equation (2.181 does not admit solutions in either L 2 
or H 1 , due to oscillations of amplitude r" 1 outside the domain of uniform ellipticity of the linear part. The 
argument is due to Johnson and Pan [To]. We truncate a solution of (2.181 at an arbitrary point, chosen to 
allow close approximation to the vortex profile Q^ m K Define, 



,-rn p(m) 



(r))=0. 



(2.19) R b = j-^2 + 2Vl + b 2 m 2 > —. 

\b\ \b\ 

Proposition 2.1 (Localized Self-Similar Profiles). Let a > Cr\ > where C > is a fixed constant and a,rj 
are sufficiently small parameters. Then for \b\ > sufficiently small, there exists G _ff 1 (]R 2 ), supported 
on \y | < (1 — rj) Rb, with the following properties. 
• Simple Profile: 



(2.20) 



(2.21) 



Q 



(m) _ 



e im0 e- ibl ^Pl m) (\y\), for P b {m > real-valued, non-negative 



(m) 



Algebraic Proximity to Q^-" 1 ^ . 

AQ< m) 



ibAQ { b m) 



Qi m) IOi m) l 2 



-*6, 



for an error term ^b, supported on {1 — rj) Rb < \y\ < (1 — rj)Rb, that satisfies the estimate, 

C(P) 

||P(y)V ^bUz, 00 < e 16 j for k — 0, 1 and any polynomial P. 
Uniform Proximity to Q( m ) : 



(2.22) 



(2.23) 



(2.24) 



y\ (q 



(m) 
b 



Q (m) )llc3 + ||e 



C\y\ 



db Qb +l 4 y 



Supercritical Mass and Degenerate Energy: 
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d{b 2 Y 



II 2 



b 2 =0 
IS 



M 2 |q 



(m)|2 



denoted d r , 



as 



and, 



0. 



< e 



-(l+Cr))(l-a) T 



The proof of Proposition 2. 1 is similar to that given by Merle and Raphael 18,19, 22 in the case of to = 0. 
An overview of the proof, and description of the particular adaptations for to ^ 0, is given in Appendix [A} 

Later in the argument, Section [2. 3[ we will introduce the linear radiation induced by the truncation error 
^b- A quantity Tb, related to the decay of this radiation, will be an important dynamical quantity, measuring 
the rate of mass ejection from the singular region. At the time we formally define r&, Proposition |2.3[ we 
will also prove the following estimate, 

(2.25) e -(i+cn)f < r b < e -^- Cri ^. 
2.1. Decomposition &2 Modulation. 

Lemma 2.1 (Modulation Near Q^). Suppose that v G H}* is close to Q^ m \ up to symmetries: 

(2.26) ^ = ift: ,) + ^ »)(^ e " ,7 ^ 

for some symmetry parameters X v > 0, b v > and j v G K such that the error is comparably small, 



(2.27) 



/ 



V y £v(y)\ 2 dy + 



-\y\ 



dy<T 



\y\<-, 

7 



6„' 



where y denotes -f- , and such that the deformed profile is sufficiently close to Q^ m \ 



(2.28) 



A„ < Yq^ v an ^' ^ v ^ a * 



Then there are parameters Xq > 7 60 > and 70 G R, nearby in the sense, 

,A 



r -i|<r| , 



(2.29) |fe - & 
and smc/i that the error eo corresponding to these parameters 

(2.30) eo(l/) = Ao«(A y)e*» 
satisfies the following orthogonality condition^ 
(2-31) Rc(e 0j |y| 2 Q^ 



2 n(-)\ l m / eo , A 2 Q^ ) = Im ( e , AQ^ 



= 0. 



Let us reiterate and extend the notation alluded to by equation (1.13) 

' e) l 



t\—e He I e 



e 2 = e im9 Im (e~" n9 - 



(2.32) 







Re I e- lb ^P, 



5(m) 



im8 



Im e 



(m) 



i6. 



Products between the components of e and behave as if they were real- valued, as does the modulus, for 
example, |e| 2 = |ei| 2 + |e2 1 2 - Moreover, since \y\ 2 and the scaling operator, A = 1 + y- V B , are radial operators, 
the algebraic relations for \y\ 2 Q^ and AQ^ are exactly the same as the case m = 0, [19, Proposition 9 

, is true regardless of m. This 



4 m) (A ei )-A(i^ ei ) 



(m) 



(iii)]. In particular, one may verify that, (ei) 

is the essential relationship for Lemma 2.2 below. In the notation of (2.321, the orthogonality conditions of 
equation (2.311 can be written, 



<e 1 ,|y| 2 I]) + <e 2 ,|y| 2 e) =0, 
<e 2 ,A 2 S)-(e 1 ,A 2 6) =0, 
(e 2 ,AE)-(ei,Ae) =0. 



These are exactly the same form as in the case m = 0. Indeed, the proof of Lemma 2.1 an implicit function 
argument, is identical. See 28 Lemma 2] for a clear exposition. For m — 0, the following Lemma was proven 
by Merle and Raphael [l9j equation (116)], and the same proof applies here. 

Lemma 2.2. Let e £ H} m )> an d assume the Spectral Property is true. Then, 



(2.33) 



^1 eijEi 



;i ,LV m) A 2 Q(™))( ei ,AQ(™)) 



|AQ(™)|| 2 2 



> 



IV„d 



|e| 2 e-l y l 



ei,<5 



(m) 



ei,|»rQ 



2 n (m) 



Definition 2.2 (Description of Initial Data). Define V^ m ) to be those functions uq £ H} m ) f or w hich there 
are parameters Ao > 0, 60 > an d 7o € R satisfy the following conditions. Let eo denote the error in 
approximating u with these particular parameters, 



(2.34) 



We require that the orthogonality conditions (2.31) are satisfied, that there is, 



^ These orthogonality conditions were introduced |l9[ Lemma 6] , and lead to a better estimate on the phase parameter than 
achieved in 181. 



proximity to Q^ m \: 



(2.35) 



in L 



Q<bl + \\e\\l 2 <(a*) 2 , 
V y e (y)\ 2 dy- 



\y\<> 



parameters consistent with the log-log rate,: 



(2.36) 



normalized energy,: 



(2.37) 



< An < e 



A^ |<ri o °. 



ko(y)| 2 



and, 



-\ y U y < r 



b ' 



Remark 2.3 (V (m > is Non-Empty). Choose b and A to satisfy (2.35) and (2.36). Let / e ffK satisfy 
orthogonality conditions \2.Zl\ with \\f\\ H i =1, (f,Q {m) ) = 1. Such an / may be computed explicitly from 



g(" 1 ). Note that d v E[Q^+vf\ 



the order of E[Q ( b m) ] to satisfy (|2.24|). 



(JF, Q( m ^) — — 1, and therefore we may choose eg = i// with j/ of 



For t he r emainder of this sect ion, we consider a fixed representative uq € r p( m \ By the continuity of the 
flow of (1.1) in H 1 , and Lemma 



2.1 



there exists continuous functions A(i) > 0, b(t) > and 7(4) € K and 
some maximal Thyp € (0,T max ] such that the following relaxations of (2.35), (2.36) and (2.37) hold for all 
te[0,T hyp ): 



(2.38) 
(2.39) 

(2.40) 
(2.41) 



\V y e(t,y)\ 2 dy 



0<b*(t) + \\e(t)\\l 2 <(a*)r, 
\<t,y) 

■ 10 

-Wi 

7T 1 

b(t) . . _ e 10 b(t) 



dy<r 



6(t)' 



< A(t) < e" 



and, 



A 2 (t)|£;o| < rg {t) . 



2.1 



at any t S [0,Th yp ). Therefore, 



Note that as a consequence of these hypotheses, we may apply Lemma 
one of the following occurs: 

Case 1:: Th yp < T max , and one of the hypotheses fails at t = Thyp, or 

Case 2:: Th yp = Tmax, 6 — > as t — > T max , and due to (2.40) we have blowup. 

In this section we will show that Case 1 cannot occur. Then, assuming Case 2, we will derive the conclusions 
of Theorem 11.11 

Remark 2.4 (Parameters). The parameter rj > 0, already introduced, rela tes to the cutoff and shape of the 
singular profile Q[" 



Parameter a > 0, to be introduced in Section 



2.3 



will be related to a cutoff point 



of the linear radiation associated with Q b . The value of 77 is determined by the value of a so that the 
argument of Subsection |2.4.2| is successful. At all times, a* > is assumed sufficiently small for all the 
appropriate constants to cooperate. 

2.2. Conservation Laws &: Basic Estimates. By substitution of the time-dependent version of the geo- 
metric decomposition (2.34), the conservation laws of (1.1) and the orthogonality conditions (2.31) lead to 
some basic estimates. 



Lemma 2.5. For all t e [0,T hyp ), 
due to conservation of mass, : 

(2.42) 



b 2 + / |e| 2 <KT 



due to conservation of energy, : 



2Re <Je, Q ( ™> - ibAQ 



(m) 



|V y e| 2 dy-3 



|Q (m) £i| 2 



(2.43) 



IQ (m) e 2 | 2 , 



Proof. To prove (2.42), expand the conservation of mass, 



IQi m) | 2 ^~ / |Q (m) | 2 ^ + 2i?e( £ ,Qi m) )+ / |e| 2 = / K| 2 - / |Q (m) | 2 dj/ 



Recognize ||jr2, and recall (2.23). Use initial condition (2.35) and hypothesis (2.39). To prove 

(2.43), expand the conservation of energy, as in 19 eqn (188)]. Use the normalized energy (2.41) to estimate 
X 2 Eq- For the terms C(e 3 ), use the exponential decay of Q^ m \ the Hardy-type inequalities below, and 
hypothesis ( |2.39| . □ 



Lemma 2.6 (Hardy- type Inequalities). For any k > and for all v E H 1 (M 2 

(2.44) / \v{y)\ 2 e-^\ < C{k) ( f \Vv(y)\ 2 + [ \v(y)\ 

J v m 2 \J J\y\<l 



2 e -\v\ 



(2.45) 



\v(y)\ 2 <Ck 2 log k[ \Vv{y)\ 2 + / \v(y)\ 2 e'^ 
|y|<« \J J\y\<i 



Proof. Equation (2.45) is proven 22 equation (4.11)], and the same techniques prove (2.44) 



Let us reiterate the notation y = £ R 2 , and introduce a rescaled time, 



(2.46) 



s(t) 



4 d- 



o A 2 (r) 



so where sq — e 4b o and, si — s(Th yp ) <E (so,oo] 



In these new variables, equation (1.1 ) now reads 



JLq(»0 + lCs - M (e) + l bAe =i(^+b) AQ^ + ^Q { ™ ] 



(2.47) 



/ ( y + b \ Ac + 7 s e 
9 b -R[e], 



□ 



where we introduced the new variable, 'y(s) = — s— 7(s), the term R[e] corresponds to those terms of u|u| 2 that 



are formally C(e 2 ), and M is the linearized operator near Q ™ , analog ous to LM ; Hl.iob. Using our choice 



of notation ( 2.32[ ), equation (2.47) has exactly the same form as that given by Merle & Raphael p9l Lemma 7] 
in the case of m = 0. Indeed, the algebraic structure in H( m } is the same, and the arguments of |19| Appendix 



C] (or 27 Appendix A]) prove the following Lemma, without modification. 
Lemma 2.7. For all s <E [sq, s\), 



due to orthogonality with |y| 2 <3j, m '\ AQ^ , and estimate (2-4S), 



{2.4 



b\ + \b s \< T 1 -^ + ( / |V,£| 2 dy + [ \e\ 2 e~\y\ dy 



due to orthogonality with A 2 Q^ n \: 



(2-49) |7. 



;i ,iV m) A 2 Q (m) 
l|AQ(-)||l 2 



\<Tl- Cr > + S(a*)(j \V y e\ 2 dy + J^Je\ 2 e-\y\dy 



in 



In order to show the coercive control (2.391 does not fail, we prove the following Local Virial Identity. 
This estimate was originally shown by Merle & Raphael in (2l] and was inspired by the work of Martel & 
Merle 13 in a proof of soliton stability for the generalized Kortcweg-dc Vries equation. 



Lemma 2.8 (Local Virial Identity). For all s £ [sq, si), 



(2.50) 



b s > Si 




klV 



\y\ 



dy 



\y\<f 



where 5% > is a universal constant. 



Proof Outline. We begin the same as the proof of (2.48): take the real part of the inner product of (2.47) 



by AQ b ' and use (2.43) to eliminate the terms 0(e), as in 

% Ui,AQ 



19 



Appendix C]. The interim result is, 



(2.51) 



bs\\\yQ {m} \\h>n^(e,< 



(m) 



il-C?7 



5(a*) 



( j |V,e| 2 dy + / 

V J\v\<^ 



\ 2 P -\v\ 



dy 



where we have used the preliminary estimate (2.48). We have also used the proximity of to Q( 



(2.22) , to isolate the ^-dependence from interactions of the form yQ^ ') e 2 as a lower-order potential, the 

same form as 19 equation (215)], here included as part of the final term. To prove the Local Virial Identity, 
use the preliminary estimate on j s and the Spectral Property, Proposition [TTTJ as adapted by Lemma |2.2| □ 

2.3. Lyapounov Functional. We cannot hope to prove b 2 is monotonically decreasing, since it is a mod- 
ulation parameter, and thus cannot hope to control e by the local virial identity at all times. In this section 
we prove a Lyapounov functional based on the mass ejection from the singular region, to which b 2 is related, 

(2.23) , and which we expect b 2 to track. To do this, we will further approximate the central profile by 
including a linear radiative tail. 

Proposition 2.3 (Linear Radiation). Fori] > sufficiently small, and all \b\ > sufficiently small depending 
on r), there exists a unique solution C,^ 6 H^ m ^(M. 2 ) to 



(2.52) 



AC (-)_ C M 



i&AC 6 (m) =* 6 , 



where is the truncation error given by (2.21 ). Radiation Q 
exists. We denote this decay rate as, Ti,. 
• Size in H 1 and Derivative by b: 



(m) 



L 2 (R 2 ), and, moreover, lim\ y \^ +00 \y\\(^ m) (y)\ 



(2.53) ||CrilAi<r 6 1 -° , J •««! 

• Decay past the support of "Iv 

( 2 -54) \\\y\\ci m) \ + \m ,v S6 

• Stronger decay far past the support of 



,W»0| 



I^C* m) l|ci<rj 



(\y\>R b ) < 1 b < 



-Cr, 



-CO. 



e -(i+cu)f < - rb < 
5 



|y| 2 ld m) l 2 l 



L°°{\y\>R 2 b ) 



< e 



which we have already discussed, equation (2.25), and, 



(2.55) 



M 2 |vcri 



\L°°(\y\>Rl) 



r 2 

< C -±- 
~ \b\ 



The proof of Proposition 2.3 is given due to Merle and Raphael (l9j Appendix E] and (22l Appendix A] 
Brief discussion of the necessary adapatations will be given at the end of Appendix |Aj 
We denote. 



(2.56) 



A(t) = e a W). 



so that, 
li 



r; 



< a<t; 



where a > is a universal parameter. Let 4>A denote a smooth cutoff function of the region, l{| y |>2A}- The 
truncated radiation, Q 171 * 1 = (1 — <^)Q , is algebraically close to and satisfies, 



(2.57) 



AC, 



C { b m) + ibA( ( b m) 



F, where, 



\F\ L ~ + \yVF\ L ~ < -j- 



We will now repeat the calculation of the local virial identity, this time including the linear radiation C as 
part of the central profile. That is we write, 

(2.58) e = e-Cl m) ' 



X(t) 



oi m) 



^(m) 



-i 7 (t) 



without affecting the parameters. This leads to a refined version of equation (| 2.47 [) for e. The proof of the 
following three Lemmas is virtually identica^o that of Merle and Raphael, |22[ Chapter 4]. 

Lemma 2.9 (Radiative Virial Identity). For all s£ [sq, Si), 



(2.59) 



dji >S 2 ( J |V y ?1 2 dy + J^ <w l^ Hyl dy 



2 J A<\y\<2A 



\e\ 2 dy, 



where 82 > is a universal constants and, 

lM||2 



A(s) = ~\\yQr\\h + ~tel /yVQi m) C 



Mm) 



Im(e,AC, 



(m) 



In the light of estimates such as (2.45 1 we cannot expect the radiative virial identity to give a good control 

1 



for e. Let (p^ denote a smooth cutoff function of the region l{|y|>3A} with steady derivative 4>'oo ~ 3 , 
the region A < \y\ < 2A. 



on 



Lemma 2.10 (Mass- Ejection). 



(2.60) 



Wdy\> 



400 



A<\y\<2A 



e\ 2 dy-T* / \V y e\ 2 dy-T 2 



Remark 2.11. As a heuristic, assume that e « Cfc on the region, A < \y\ < 2 A. Use the definition of T b to 
approximate the mass. Then with hypothesis (2.391, Lemma 2.10 suggests continuous ejection of mass from 
the region \y\ < ^, regardless of whether that region is growing or contracting. 



Together with the conservation of mass, Lemma 2.9 and Lemma 2.10 prove the following Lemma. The 



argument relies on (2.53) and the relation between parameters a and 77 stipulated by Proposition 2.1 



Lemma 2.12 (Lyapounov Functional). For all s £ [so,s\), 
(2.61) 



d s J <~Cblr b + [ | V^ 2 dy+ [ m 2 e-M dy + [ 



A<\y\<2A 



where C > is a universal constant, 



(2.62) 



+ 2 (ei, E) + 2 (e 2 , 9) + / (1 - \e\ 2 dy 



S 2 
800 



bfi(b) - f h(v) dv + blm (e, A^ m) ) J , 



and f\ is the principal part of fx , 

Mb) 



\vQ 



(m),|2 



\L 2 



-Im 



^Where Merle and Raphael write, f;, = £ re + iCimi one should instead read, = Cl + *C2i each component with spin m 

following the convention of equation (2.32 1. 



2.3.1. Estimates on J . To first order, J quantifies the excess mass remaining in the singular region. After 
explicitly accounting for this mass, J is comparable to ||e||^i, up to a power of Tf, that depends on our 
choice of truncation of the radiation. 

Lemma 2.13. For all s € [so,si) we have the crude estimate, 
(2.63) \J-d m b 2 \ <5 3 b 2 , 

where < S 3 <C 1 is a universal constant, and d m b 2 is the approximate excess mass of profile Q{ . 
Lemma 2.14. Let f 2 denote those terms of J that are formally 0(b 2 ), 



f2(b) = \\Q i r ) \\h-\\Q (m> \\h-^ [bh(b) 



fi(v) dv 



These are the terms concerned with the excess mass. For all s € [sq, si) we have the refined estimate, 



(2.64) 



J{s)-h(b(s)){ 



< Tl- Ca +CA 2 log A I [ |V y e| 2 dy + [ 

-n-°' a +^( f\Vye\ 2 dy+ [ 

C \J J\v\<^ 



\e\ 2 e-\y\. 



e\ 2 e-^ 



dy 



Proof. The crude estimate (2.63) can be either proven directly or seen as a special case of (2.64) and 
hypothesis (2.39). The estimate (2.64) and its proof is exactly as given by Merle & Raphael, 22, equation 
(5.6)]. The essential point is that the most difficult term of J(s) — f 2 (b(s)) can be handled with the 
conservation of energy (1.3), here written in rescaled variables, 

2<ei,E)+2(e a ,e>+ / (1 - M\e\ 2 = (l^e^) + (xL m) e 2 , e 2 ) - J A |e| 2 

+ £[g[ m) ]- X 2 E + O(e :i ). 



To establish the lower bound of (2.64) we need to be coercive. We claim that, 



Lemma 2.15. For v = v\ + iv 2 € H} ■,, 



(2.65) 



4 TO V,«i) + (L ( " l) v 2 ,v 2 ) > S 3 \\vf H1 - i ( [r,.,, 



vi,VQ 



(to) 



V2,Q 



(m) 



where 0+ is the normalized eigenvector corresponding to the smallest eigenvalue of L 



(to) 



Merle and Raphael 22 Appendix D] remark that (f> + lies in the span of Q(") and \y\ 2 Q {rn \ and that 

<Pa\v\ 2 



(2.65) may be localized to, 



r ( m ) \ i / t ( m ) 



> S 2 \\v\\ 2 Hl - ~ (( Vl ,g (m) ) 2 + \y\ 2 Q {m) ) 2 + (« 2 ,Q (m) ) 2 ) , 
assuming A is sufficiently large relative to the exponential decay of Q^ m \ 



□ 



Proof of Lemma \2.15\ Following the variational characterization of Q^ m \ Weinstein 32 Prop 2.7] argues 
(in the case m = 0) that for all / G Hly d 2 J[Q^ + e/]| £=Q > 0. By an explicit calculation we conclude, 



inf (4 ro) /,/)>0. 

Let /i + < be the lowest eigenvalue of L^, and <f>+ € L 2 the corresponding normalized eigenvector. If 
there were two linearizely independent negative directions, then there would be one perpendicular to R^ m \ 
Therefore, 

inf (4 m) /,/)>0. 
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The following argument due to [14] is an improvement on the proof of [32] Prop 2.9]. Consider, 
5+=inf{(i ( f " l) /,/) I ll/IUj m) =l and (/,</>+> =o} >0. 

Assume <5+ = 0. Then by weak convergence of a minimizing sequence there exists a minimizer /+ € ^} m \, 
and there are lagrange multipliers so that, 



An inner prod uct with / + implies t\ = 0, and then an inner product by implies t% — 0. As we remarked 
Warchall 26 found that the nullspace of restricted to H} m \ is empty, and we 

□ 



Secti 



1.2 



Pego 



have a contradiction. 

2.4. Description of the Blowup. Let us consider the hypotheses of Section [2~T] in turn. In each case, we 
will show that if the solution exists for t = T} lyp , then the hypothesis holds for some interval [0, Th yp + S). 
This will prove that Case 1, introduced in Section 2.1 cannot occur, and that, therefore, Ih yp = T max . This 
means that the dynamics described by ( |2.38[ ), ( |2.39 1, (2.40) and (2.411 persist for the remaining lifetime of 
the solution. Indeed, we will show that that lifetime is finite, equation (2.181, and use these dynamics to 
prove the behaviour claimed for Theorem |1.1| 

2.4.1. Hypothesis (2.38). Preliminary estimate ( |2.42 ) shows that hypothesis (2.38) cannot fail. 

2.4.2. Hypothesis KM). 



Lemma 2.16. For all s G [so,Si), 
(2.66) 



\V y e\ 2 dy 



\e\ 2 e 



-\y\ 



dy<H 



\y\< ] 



b • 



which shows that hypothesis {2.39) cannot fail. 
Proof. Consider arbitrary fixed s € [sq,Si). 



(a) If d s b(s) < 0, then (2.66) follows from the local virial identity, (2.50) 

(b) If d s b(s) > 0, then there exists a largest interval (s+, s), with s < s. 

(c) s + 

This implies, b(s + ) < b(s) and either, or, 

(d) d s b(s+) = 0. 



on which d s b > 0. 

so, 



In case (c) or (d), either by the initial condition (2.35) or the local virial identity, respectively, 

|e( S+ ,j,)| 2 e ^l dy<Tl 



\V y e(s + ,y)\'dy + 



\y\<j 



- b(s 



From the upper bound of refined estimate ( |2.64 1, and assuming a > is sufficiently small, 



(2.67) 



J(s + )-f 2 (b(s+))<n 



< ri 



(2.68) 



Since J is non-increasing, and from the lower bound of refined estimate (2.64), 

rf w >J(s) - f 2 (b( s+ )) 

\e(s,y)\ 2 e-\y\ dy) 




l^ a + (f 2 (b(s))-f 2 (b( S+ ))). 



As noted in the proof of crude estimate ( 2.63[ ), we may assume the constant S 2 of equation ( |2.59 1 is suf- 



ficiently small relative to do, such that < 



b 2 =0 



< oo, and proving that (/2(&(s)) — /2(&(s+))) > 



0. Assuming a > is sufficiently small, this proves (2.66). 



□ 



2.4.3. Hypothesis {2.4ty . 
Lemma 2.17. For all s E [sq, si), 



(2.69) 



6(a) > 



3tt 
4 log s ' 



and, A(s) < \/Xoe 3 1o s 



Proof. Recall the bounds on T^, (2.25), hypothesis (2.39), inject into the local virial identity (2.50), carefully 
integrate, and recall the clever choice of s , (2.46), 

6, 37T 



<9,e + « = — ^— e +4 <> <1 
b 2 4 



e +4b < s _ Sq jf. e +« = s. 



Next, we take hypothesis (2.39) and preliminary estimate (2.48) to approximate the dynamics of A, 

7T 



(2.70) 



b\ + \b s \ <n 



X s 26 

— > — > . 

A - 3 ~ 2 logs 



By the initial condition on bo, (2.35), we may assume sq is sufficiently large so that, 2 lo a ^ a — 



By the initial choice of a log-log relationship, ( 2.36 ) , we may assume, — log Aq > e 2b o = s 



3 \^logs log s 

That is, by integrating (|2.70|) we have, 



, , 1 . . IT S 

logA > --logAo + . 

2 3 log s 



□ 



Corollary 2.18. 



Corollary 2.19. 



I hyp 



A 2 (cr) da < A 



e ^ da < a*. 



A < e 



, which shows that half of hypothesis (2.4ty cannot fail. 
igain assumin; 
log(sA(s)) > 



Proof. Due to equation (2.69), and again assuming so > is sufficiently large 

7T S 



3 log s 



log s > 



logs' 



Take the logarithm and apply equation ( 2.69[ ) again 
(2.71) log|-log(sA(s))| >lo. 



logs 



4 , 7T 

> — log s > 

- 15 & ~ 56(a) 



sA(s) < e 



-e56 



□ 



Lemma 2.20. For all s <G [sq, s\), 
(2.72) 



6(a) < 



An 
3 logs 



Proof of Lemma \2. 20[ Due to the crude estimate (2.63) and the Lyapounov inequality (2.61 1, 

9,e+ 



d ^ b 5tt , / d 

J > jF b e 4 V a > 1. 



where the final inequality is due to | > 1 + Cry, the bound for Tb, (2.25), and assumes a* is sufficiently 
small. By integrating the inequality, 



(2.73) 



e T 4 v J(s> > e 



5tt / ^0 
' 4 V -7< s o) 



S - Sq. 



Finally, by the crude estimate (2.63) and the definition of s (2.46), 



57T / d TT 
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which, again with the crude estimate ( 2.63 ), proves ( 2.72 1 from (2.73). Finall y, we note here a related 
estimate that will be used in Subsection 



2.4.5 



Divide the Lyapounov inequality (2.61) by V J ', integrate in 



time, and use the crude estimate once again to get, 



(2.74) 



Corollary 2.21. 



r 6(CT ) + / |V,e| 2 dy + ^ iq |e| 2 e-^ dyj da < ^J{^) - VW) < bo- 



< A, which shows the other half of hypothesis (2.4-0) cannot fail. 



Proof. From the approximate dynamics of A, equation (2.70), 
A 



47T 

-f < 36 < 

A logs 



f 1 

— log A(s) < — log A + 47r / da 

J s (> logtJ 



Bound the integral with 4ir(s — so), and by using (2.72 ) again, e 47r ' s s °- ) < e 
of so (2.46) and initial condition (2.36), 

A(s) > A e 4 ™ e -^^ > 



4tt( e 36 ( s ) —s 



□ 



With the definition 



□ 



2.4.4. Hypothesis (2.41). As another consequence of the approximate dynamics of A, equation (2.70), 

'A, 



d_ 

ds 



< - Ybe^b < 0, 



A 2 (t)e"< t ) < Ape^o . 



Then, due to the bounds on Tb, \2.25 1, and the initial condition ( 2.37| ), 

x 2 (t)\E \ < r 4 (t) e^x 2 \E \ < r 4 (t) e^r™ « r 4 (t) . 

This shows that hypothesis ( |2.41| cannot fail. 
2.4.5. Dynamics of Theorem ] 

Proof of Log-log Rate. By proving Thyp = r max , we have already shown blowup in finite time, due to Corol- 
lary [2T8J Here we establish the rate. By direct calculation and a change of variable, 



(2.75) 



-d t (A 2 log|logA|) 



■log log A 2 



|logA|log|logA| 



Recall the approximate dynamic -jRid, and with hypothesis (2.40), equation (2.75) reads, 

~ <-9i(A 2 log|logA|) <C. 
Integrate over [t,T max ). Since A is very small we may estimate, 



(2.76) 



1 



T — t 

1 max v 



C Vlog|log(T max -t)| 



< A(t) < C 



log|log(T max -t)| 



Moreover, the relationship between X(t) and the log-log rate has a universal asymptotic value as t — > T max , 
see [22l Prop 6]. □ 
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Proof of Singularity Description in L? . The proof here is heavily inspired by that given by Merle and Raphael, 
[20| Section 4]. First, we show for any R > there exists u* such that, 

(2.77) u(t)^u* in L 2 x (\x\>R) as t -> T max . 

Second, to establish equation ( 1.12| ), we will prove that both, 



(2.78) 



u* £ L 2 (R 2 ), and, 



= lim 

t->-T m . 



K*)l s 



Let e > be arbitrary. Choose some T ma x — e o < t(R) < r max . By hypothesis ( 2.40[ ) we may assume that, 
it(t) = u on {|x| > ^} for t 6 [t(-R), r max ) and by equation (2.74) we may assume that, / t /m* J |Vu| 2 dxdt < 
eo- For a parameter r > 0, to be fixed later, we denote, 

(2.79) v T (t,x) =u(t + T,x)-u(t,x). 

Since t(R) < T max , u(t) is strongly continuous in I? at time t(R). Thus, there exists r such that, 



(2.80) 



J \v T (t(R))\ 2 dx < e for all r € [0,r ]. 



Denote 0# a smooth cutoff of the region Inj^fl}. By direct calculation, 
1 

(2.81) 



9t / <Pr\v t \ 2 =Im 



Im 



V ' 4>r ■ V/v r dx 

4> R v T (u\u\ 2 (t + t) - u\u\ 2 (t)j dx 



Regarding the first RH term of (2.81), use Holder, (2.80), and the choice of t(R), 

i 



t(R) 



Im / V0 fl • Vv T v T dx ) dt\ < C 



t(R) 



V dt e§ < Ce . 



|fi|li« < Nl? ' — 

to get control by e . We have proven that u is Cauchy on \x\ > R, 



Regarding the second RH term of (2.81), control with \\u\\^4, < \\u\\i J 2\\u\\ 2 ^ 1 <C and integrate in time 



J 4>R\v T {t)\ 2 dx < Ce for all r G [0, r ] and t € [t{R),T„ 



r). 



We now turn our atten tion t o (2.78). The profile and radiation have support of radius R(t) = A(t)X(t), 
which, due to hypoth esis ( 2.40), is going to zero with a bound, R(t) < (T max — t)z~ s . From the definition 
of A(t) and equation (2.45) we may bound J(l — 4> R(t))\u{t)\ 2 to prove, 



<f>R(t)\u(t)\ 2 -- 

which proves that the following limit exists, 



lim 



lim 



lim 

*->-T ma3 



4>R{t)\u{t)\ 2 = 

d?R(t)Mt)\ 2 



lim 



|5(*)| J 



This completes the proof of (2.78) and (1.12). 



□ 



3. The Spectral Property 



We now provide a numerically assisted proof of the spectral property for the case m — 1. We also present 
some computations on higher order vortices and discuss why they do not work. Before proving the Spectral 
Property of Section [L4} we will establish the following variant: 



Proposition 3.1. Let e € H} m \ satisfy the orthogonality conditions, 



(3.82) 



ei,<2 



(rn) 

M\ - 



(m)\ _ 



C-2 



AQ< m >) = (e 2 A 2 Q [m) ) =0. 



Then, for the case m = 1, £/iere is a universal constant C m > 0, so £/ia£. 



(3.83) 



ft 



(m). 



(e,e)>C m / (|V„ e | 2 + e-l»l|e| 2 )dj/ 
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Proposition 1 1 . 1 1 is an immediate corollarjJ^J Following [6|[l5], we proceed in two steps. First we count 
the number of negative eigenvalues of the operators and C^ 1 . We then show that the assumed L 2 

orthogonality conditions are sufficient to project away from the negative directions of the bilinear forms, 
and H^, associated with and C2 ■ 

We now restrict ourselves to e £ H(m)> e = e "™ e rad ; where e ra d is a purely radial function, 

(3.84) e rad G H^ d+ = Hl d (R 2 ) n {u \ {x^u G L 2 (R 2 )} . 
Given e € H}-. , we calculate 

(3.85) Ci°»e = e«* (-£ - \± + £ + SR^y ■ V«±) e rad 
This motivates defining the two operators and inner products on H^ &d _^_ 

2 2 

(3.86a) = -A rad + ^ + 3R^y ■ VR {m) = -A rad + ^ + Vi, rad 

2 2 

(3.86b) £g d = -A rad + ^ + • Vi?( m > = -A rad + ^ + V 2 , rad 

(3.86c) <> d (, •) ee (£g d , •) , < ad (-, •) - (42d«. ) , 

where A rad is the radial Laplacia n, A r ad ee £^ + 

The orthogonality conditions, (3.82 1, are now formulated as 

(3.87) (efSitf™^ = /ef d ,Ai^ m )\ = <Je^ d ,AE( ro A = <^ ad , A 2 R^ m A = 0, 

where 

e = e lme (e^ + ie^ ad ) 

All that follows relies on the reduction to a series of one dimensional radial problems. 
3.1. The Index of Bilinear Forms. 

Definition 3.2. The index of a bilinear form B with respect to vector space V , denoted in&yB, is the minimal 
co-dimension over all subspaces of V on which B is a positive. 

For bilinear forms induced by self-adjoint operators {i.e. B = (£■,■)), the index corresponds to the 
number of negative eigenvalues of the operator. To calculate the index, we extend Theorem XIII. 8 of Reed 
& Simon f29l to: 



Theorem 3.1. Let U solve, 



d 2 1 d m 2 
CU = --^U--^-U + V(r)U+ r ^U = 0, 
dr A r dr r z 



with initial conditions given by the limits, 

lim r- m U(r) = 1, lim -f- (r^ m U(r)) = 0, 
r-yO r-yO dr y ' 

and where the potential V is sufficiently smooth and decaying at oo. Then, the number of roots of U not at 
the origin, N(U), is finite and equal to the index of the bilinear form associated to C over the vector space 

TTl 

n rad+- 

Proof. The proof, which we omit, is quite similar to the proof of the indicated Theorem of Reed & Simon. 
In turn, that proof is a generalization of the Sturm Oscillation theorem for two point boundary value 
problems. □ 



! See the end of Subsection 3.3 for details. 
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Figure 3. Plots of the functions and Z (m) solving ( |3.89[ ) for m = 1, 2, 3. The number 
of non-zero zero crossings determines the indexes of H\ and H^- 



Proposition 3.3 (Numerically Verified). For the cases m = 1,2,3, 
(3.88) 

Proof. Using the methods described in Appendix 151 we solve 



£\lLu (m) = 0, lim r- m U {m) (0) = 0, lim — r- m U {m \r) = 0, 

i ' Iad r->-0 V ' r->0 dr 

£ { 2 ™i d Z^ = 0, lim r - m Z^ (0) = 0, lim —,•-"»#("») (r) = 0. 



(3.89a) 

(3.89b) „ .... . . 

Plotting the solutions in Figure [3} we can see that I7( m ) has two zero crossings and has one zero 

crossing. Subject to the acceptance of these computations, Theorem |3 . 1 1 yields the result. □ 



Proposition 3.4. There exists a constant, 5q > 0, depending on m, such that for 5 € (0, 5q), the bilinear 
forms associated with the perturbed operators, 



C 



(m) 
j.r&d 



c 



(m) 
j.rad 



Se 



-\v\ 



have the same index, i.e. 



(m) 
j,rad- 



Proof. We briefly sketch the proof, which follows from three observations. First, the of the solutions of the 



perturbed form of (3.89) are continuous with respect to 8. In particular, there is Cy oc convergence. Second, 
the roots of the index functions, in the perturbed and unperturbed cases, must be simple. For a sufficiently 
small 5q, we can ensure that on any compact interval the perturbed and unperturbed solutions have the 
same number of zeros. Finally, for a sufficiently large compact interval, outside the interval the equation 
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Table 1. Inner products associated with L\ r£ [ d for different winding numbers. 



in 




K { 2 m) 


^3 


K (m) K (m) _ (gWy 


1 


-0.48237 


-25.798 


1.28129 


10.8025 


2 


0.520152 


-13.1545 


1.7983 


-10.0762 


3 


2.59249 


5.1232 


-1.54694 


10.8888 



Table 2. Inner products associated with C 2 r£ [ d for different winding numbers. 



m 


T (m) 
J l 


T (m) 
J 2 


T (m) 
J 3 


T (m) j(m) 


-{4 m) T 


1 


6.6985 


163.548 


-47.7764 


-1.1871e-| 


-03 


2 


25.1685 


1319.28 


-235.186 


-2.2108e^ 


-04 


3 


82.6396 


8426.22 


-936.752 


-1.8116e4 


-05 



is approximately "free" (the localized potentials are negligible), and we can ensure there are no additional 
zeros; this may require further shrinking Sq. □ 



3.2. Orthogonality Conditions and Inner Products. To verify that orthogonality conditions (3.87) 
project away from the negative subspaces, we need to compute a number of inner products of the form 
jTad"'")' where u solves C^^u = /. Although these products are computed numerically, we justify 
their existence: 

Proposition 3.5 (Numerically Verified). Let f be a continuous, radially symmetric, localized function satis- 
fying the bound \f(r)\ < Ce~ Kr for some positive constants C and k. There exists a unique radially symmetric 
solution, 

£j,rnd u = fi 3 = Ij 2- 

that belongs to the class, u £ L°°([0, oo)) n C 2 ([0, oo)). 



Proof. This is Proposition 2 and 4 of 6 , along with our computations of the indexes in Lemma 3.89 Sec 15 



for some additional details and a full proof in dimension d = 1. □ 



Remark 3.2. The solutions in Proposition |3 . 5 1 may not vanish as r — > oo. Indeed, they can only be expected 
to be bounded. 

Proposition 3.6 (Numerically Verified). LetU\, U 2 , Z±, and Z2 be L°° radially symmetric functions solving, 

(3.90a) 4™L c/ i = - R(w0 > 
(3.90b) = AR (m) , 

(3.90c) C™ A Z 1 =hlfr»\ 

(3.90d) £ ( ™l d Z 2 = A 2 R( m \ 
Then the inner products, 

K[ m) - (d&FuUx) , KjT> = (C^IU 2 , U 2 ) , Kt> ^ (£^Uu U 2 

T (m) _ I Am) 7 y \ r(m) _ I Am) 7 y\ j(m) _ I Am) 7 7 

J l = \^2,rad^l' Z l/ ' J 2 = \ i '2,rad Z 2, ^2 J , J3 = \^2,rad Z l > Z 2 

take the values given in Tables^ and^ 

Proof. Using the methods described in Appendix [B] these are computed numerically. □ 

As with the indices, we have stability of the inner products with respect to perturbation by a small 
portential: 
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Proposition 3.7. Let U i and Zi denote the solutions and k\ ^ and j\ "* the inner products, analogous to 
those of Proposition 



3.6 



for the boundary value problems with the perturbed operators, £ ?rad . For S > 



sufficiently small, the solutions and inner products are continuous with respect to 5. 
Proof. This follows from the invertiblity and continuity with respect to 5 of the operators. 



□ 



3.3. Proof of the Spectral Property. We are now able to prove Proposition 3.1 The arguement closely 
follows the proofs found in j6j[l5] . The two bilinear forms, and % 2 m \ are treated seperately. First, 

we will show that L? orthogonality to and AQ'" 1 ' suffices to project away from the negative subspace 

of T-Li- This will only be successful for m — 1. Later, we will show that L 2 orthogonality to AQ^ and 
\2Q(m) p ro j ee t; S away from the negative subspace of W% ■ 



Spectral Property for Ti^ . Given an element u G H} m y u = e lme it rad , satisfying orthogonality conditions 

(3.82), showing positivity of H^ 71 ^ on such a u is equivalent to showing posiviity of 'Hj I ™ ad on w rad G H^ ad+ 
satisfying orthogon ality con ditio ns (3.87) 
By Propositions 



3.3 



and 



3.4 



r(mT 



H l rad has a two-dimensional subspace of negative directions. Recall the 

notation of equation (3.90). Let V = span { U\ , U 2 } . We will prove that, for m = 1, T^ilid 1S negative on 
all of V. Indeed, consider an arbitrary element of this space, 

U = c 1 U 1 +c 2 U2, 

and compute, 



^l.radi^) = C l^l +2c lC2 K 3 +C 2 K 2 



(3.91) 



Ci c 2 



(m) 



K 



(m) 



If the above matrix is negative definite, then the bilinear form is negative on the two dimensional space 
V. We examine the matrix using the computations in Table [I] and elementary properties of matrices. For 
m — 1, 

tr = k[ 1] + K 2 1] = -26.2804 + o(l), 

where o(l) corresponds to taking the perturbation parameter, 5, sufficiently small. Therefore the sum of the 
two eigenvalues is negative; at least one is negative. Next, 

det = K { l ] K 2 1] - (K^) 2 = 10.8025 + o(l), 

so the two eigenvalues have the same sign. Therefore H^l^ is negative on V. Table [l] shows that this is 
false for m — 2,3. We restrict our attention to m = 1. 

I 2 ), we could decompose the space as 



Pretending that V C -ff rad+ ( 
(3.92) 



u 1 



(R 2 )=V( S (1) V 1 - 



(1) 



where our notation indicates that we have formed the orthogonal complement with respect to the H 1 rad 
bilinear form. The non-degeneracy of the matrix (3.91) justifies this decomposition. 

It follows that Hi Illd is positive on V- 1 . Otherwise, there would be W € V 1 - such that Ui ia d(W, W) < 0, 
which implies by construction that, span (W, U±, U 2 }, is a negative definite space of 'Hj 1 ,'' ad with dimension 
three. But then, given any subspace U C H^ad+ 01 codimension two, U fl span {W, U\, U 2 ] ^ 0, which 
contradicts the index calculation. 

Finally, given any function u € H^ ad+ and L? orthogonal to and AIv- 1 ', we decompose u as 



u = ciU i + c 2 U 2 + u 
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± 



where, u 1 - E V , again in the sense of (3.92). Then, 



= (u,RW 



L 2 



7(1) 



(1) 



ciHi; ad (f/i, u x ) + c 2 n\^ d (u 2 , ui) + Hi,; ad (u x , ui) 



L 2 



Cl 



<Vi,Ai?(" l ^ 2 +c 2 (u 2 AR [1) ) L2 + (u\Ai?W 



L- 



cin[% d (U 1 , U 2 ) + c 2 H ( li d (U 2 , U 2 ) + u^l^, U 2 ) 
CiA 3 +c 2 A 2 ■ 



Due to the non-degeneracy of (3.91 ), the only solution is ci = c 2 = 0. Therefore, for all such u, 



This yields the positivity of %^ on Hhy 

Of course, f/i and U 2 are raoi in H^ d+ . The above argument is made rigorous by introducing an ap- 
propriate cutoff function and then taking limits. We refer the reader to [6 15 ; we will not reproduce this 
here. 

□ 

Spectral Property for ■ As in the case of , we will prove positivity of H™ subj ect to th e orthogo- 



nality conditions, by working with the associated radial form, T-L 2 rad- By Propositions 

one ne 
Define 



3.3 



and 



3.4 



^2, rad naS 



one negative direction. Examing Table |2| neither Z^ nor Z 2 m ^ appears to point in the negative direction. 



(3.93a) 



(3.93b) 



Then C ( 2 m) Z = and, 



7 (m) 

R^ = AR^ - -^—A 2 R i ~ m \ 
J 2 



Z = Zi - 



J 



(m) 



(m) ' 



(3.94) 



n { 2 m \zj)= 1 



(m) 



J 



J 1 J 2 



cry 



< 0. 



Now that we have constructed a negative direction, we apply a similar argument as in the case of 'Hi™^; 
however, this will hold not just for m = 1, but also for m = 2, 3. We decompose ^ad+l^ 2 ) as 

(3.95) H^ d+ (R 2 ) = span {z} e^w span {z}^ 

Since the index of "H 2 7 a d i s one i we are assured that it is positive on span ■ Finally, given v € H^ ad+ 

orthogonal to Ai?( m > and A 2 R { - m \ it may be decomposed as v = c\Z + v , and, 

0=(v, R {m) 



L 2 



c 1 H { 2 n \z,Z)+H i 2 m \v ± ,Z) 
Cl H 2 m) (ZJ). 



Invoking (3.94), this implies that, v — v 1 - G span . Therefore, for such v, 

'^ < Cld(v,v) > 

for m = 1, 2,3. Posivitiy of H 2 n) on HL,, subject to orthgonality to AQ^ and A 2 Q i - m \ follows. 



□ 



Proof of Proposition 3.1 Given e — e± + ie 2 satisfying the orthogonality conditions of Proposition 3.1 
have proven that, 

n (1 \e,e) =H ( 1 1 \e ll e 1 ) + U ( 2 \e 2l e 2 ) > 0, 

from which we infer, 

H (1) (e,e) >S J e- M \e\ 2 dy. 

Let 9 > 0. Then, 

(l + 9)H (1) {e 1 e)>9 J\Ve\ 2 dy + 9 J Vi| ei | 2 + V 2 \e 2 \ 2 dy + 6 J e -^\e\ 2 dy. 
Although the potentials are sign indefinite, for 9 sufficiently small, 



we 



(3.96) 

We now have the result 



J V l \e l \ 2 + V 2 \e 2 \ 2 dy + 5 J e~^\e\ 2 dy > & - J e^\\e\ 2 dy. 



\Ve\ 2 dy 



2(1 



-|»IU|2 



e\'dy 



>S J |Ve| 2 + e-l^|e| 2 d2/. 



Proof of Proposition [177} Let e £ H^(M. 2 ) with e = t\ + ie 2 , and further decompose this as: 
(3.97a) ei = e t9 (u + Cl i? (1) + c 2 Ai? (1) " 



(3.97b) e 2 = e l6 (v + diAitf 1 ) + d 2 A 2 R^ 

where u _L L 2 Ai?W and w _L L 2 , A 2 i?W . Expaning, 

W< m >(e,e) = ^ m) (e 1 ,e 1 ) + ^ m) (e 2 ,e 2 ), 



M"°(<^i) 



ft 



u) + 2 C1 (£^ U , + 2c 2 (4^, AflH 
c 2 M[ m] + cjMf 1 + 2cic 2 M 3 (m) , 



□ 



H 2 m) (e 2 ,e 2 ) = n^l(v,v)+2d 1 (c^lvAR^)+2d 2 (ct\A 2 R^) 

+ d\N[ m) + d 2 A 2 (m) + 2d 1 d 2 A^ m) , 

where Mj m) ,ivj m) are fixed terms arising from applications of the bilinear forms to combinations of 

J? (m) ) Ai?M, and A 2 R ( - m \ 

We now construct a lower bound. Let 9 > 0. Then 



Cl 



(3.98) 



< 1 [9- 2 c 2 +9 2 (c ( ™i d u,R^) 
9- 2 c 2 + 9 2 (jMld&lfi 
9~ 2 c\ + 9 2 ( f | u ||£^ d i?(™)| 1 / 2 |£( 1 )i?( m )| 1 / 2 



- 2 
1 

< - 

- 2 

< C 



9- 2 c\ + 9 2 f\C^ aA R^\\u\ 
< C (V 2 c 2 + 9 2 J , 



e-M\u\ 2 



2:t 



The other terms in which u or v appears once are similarly controlled. Therefore, 

n^(e, e) > n^l(u, u) + y^Uv, v) - ce-\ c \ + cf + d\ + d 2 ) ce 2 / e -i»i | u + 

For the case m = 1, we apply Proposition |3 . 1| to get 
HW(e,e) >C {1) f \W(u + tv)\ 2 



+ (c ( i) - ce 2 ) / e -i»i |u + iu| 2 - <70~ 2 ( c 2 + c 2 2 + d\ + 4) 



> 



C(i) 



|V(lt + iv)\ l + e~\ y \ \u + iv\ l - CQ-'\6{ + ci. + dj + d 2 ), 



where we take 6 > sufficienty small. Finally, 

\Ve l9 {u + iv)\ 2 + e~M\e ie {u + iv)\ 2 



> C 



Ve 



0{c\ 



d\ 



4) 



a 



Appendix A. Almost-Self Similar Profiles 



In this Appendix, we outline the proof of Proposition 2.1 showing modifications of the proof given 



ib— r^( m ) 

e 4 Q\ ' 



Am6 p(m) 



in the case m = 0, [18 19 22 . We then briefly discuss the proof of Proposition 2.3 Recall that for 



Pi '(r) we have equation (2.181 



AP, 



(m) 



P 



(m) 



0. 



r 2 4 

This is not a scale-invariant equation, and there is no clear representative solution. Fibich and Gavish [1] 
chose to consider the solution where the boundary condition lim r _yo r^ m pjf n \r) ^= is chosen to minimize 
the amplitude of the asymptotic oscillation. Since we intend to truncate anyways, it is more convenient to 
choose boundary conditions, 

V for < r < (1 - rf)R h , 
= for r = (1 - r))Rb. 



(A.99) 



p M (r) 



Recall that was chosen, (2.19), so that the strong maximum principle applies to, A — ^1 
on a region larger than, r < (1 — rj)Rb. 



bl r 2 
4 ' 



Step 1: Existence of P, 



(m) 



Following the argument of 18 p605-606], let P( m ) denote the space of radial profiles of functions in Hh 



That is, radial H 1 functions f(x) for which x f(x) £ L 2 . Perform a constrained minimization of, 

b 2 



(m)- 



2J b [w] 



|Vw| 



\w\ 2 + m 2 



\rw\ 



over the subspace of finite-variance functions in P( m ) with u>((l — f])Rb) = and J \w\ 4 = 1, where all 
integrals are taken over a larger compact set, for example r < (1 — ri 2 )Rb. Note that J& is coercive on 



(A.100) 



Mw] > C( V ) 



IVtul- 



This minimizing sequence can be assumed to converge weakly in H}y which is simply a subspace of iP(IR 2 ), 
and thus strongly in L due to Sobolev embedding on a compact domain. Here we use that equation (1.1) is 



energy subcritical. The Lagrange multiplier of the Frechet derivative shows that (2.18) is satisfied. Interior 



regularity estimates show that the weak limit is C 3 on r < (1 — T])Rb- The weak limit is also strictly positive 
due to w((l — r))Rb) = and the maximum principle. 
Step 2: L°° Estimates, Uniform in b. 

There exists a fixed constant C > for all |6| > sufficiently small so that, 



\pr\L°°<c. 



(A.101) lp(m)| 
Moreover, there is uniform decay of the tail of the solutions. For the same b, 



(A.102) 



SU P \ p b |i°°(r>K) — > as R 
|6|~o 



Both bounds are proven in 18 p606]. Equa tion (A. 101 ) is a simple consequence of the variational characteri- 
zation of Step 1, whereas to prove equation (A.102 ), truncate to r > R, treat \pj; m ' | as a one-dimensional 
function, and control by the standard Sobolev embedding H^(M.) <-» L°°(M). 
Step 3: Local Convergence to R^ m ~> (in C 3 ). 

As b — > 0, pf- m > converges weakly to some positive radial function P, with decay to as r — > +oo, and 



23 



which satisfies, AP — 

(l + ^) P + P\P 2 \- This characterizes P as the unique groundstate i?( m ) 

Moreover, due to interior regularity estimates, on any compact set the convergence of P^ m ' is strong in C 3 , 
up to a subsequence in b. 

Step 4: Uniform Convergence to R^ (in C 3 with exponential weight) 

p658-659]. Consider the operator K, = A - ( 1 + 



Here we adapt the argument of 
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■ 2 4 ' 2 J ' 

which satisfies the maximum principle on 1 < r < (1 — ri)Rb, for any 77 > sufficiently small. Restate (2.18) 
as, 



(A.103) JCP b {rn) = |-P b (m) - (P 6 (m) ) . 

Consider the new function /(,(r) = e ^ J) - l " Rb0 (- n b) ; with, 

Note the dependence on m. By direct calculation, 



f^JCfb = (1 - 7 ? ) 



(i) 



+ (i-'?) 2 1- 



P& 



r 

P& 



- 1 



We now approximate each term on the region ^ < r < (1 — rj)Rb, 



b 2 2 V 2 
4 r T 



/ b -^/ b < ^^ + ((1-^-1) (l 



r 

P& 

2 



Recall that, Rb — V 2 + 2 \Zi+ b '» _ gy assumm g 6 > is sufficiently small with respect to 77, we conclude 
fb~ l K-fb is strictly negative for the given range of r. 

From Step 2, and the exponential decay of R^ m \ there exists a fixed value vq > ~ such that for all > 
sufficiently small, 



2 ^ 



P. 



(m) 



> for r G O = r < r < (1 - T))Rb- 
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We have shown that JC (cfb — P^ n ^ < for r 6 D, and any arbitrary constant c > 0. Now we note that, 

lim/ 6 (r ) = e- (1 -" )ro > 0, 



b-S-0 



so that we may choose our constant c = 2i?( m )(r )e+( 1 -'') ro and, with our boundary condition ( |A.99| , 
conclude that, 



ch(r)-Pi m \r) 



on 



> 0. 



The maximum principle may now be applied. The same argument can be a pplied to R^ m \ b — 0, and the 
weight f(r) — er^ r ^ r ' . With Step 3, this proves the first precursor of (2.22 1, 



(A.104) 



C 3 



as b -> 0. 



To prove the bound for the energy, (2.24), note that without loss of generality (1+Cr/)(1 — a) = (1 — S) < 1. 
Introduce a new operator /C and function fa in terms of 6 in place of and argue Step 4 again. In particular, 
we may assume that r < SRb <C (1 — ?7) 2 i?b < r < (1 — rf)Rb <C (1 — S)Rb- 
Step 5: Uniqueness of P b , Continuity in b 

For fixed 6q > sufficiently small, and b m bo, consider, 



(A.105) 



Tb,b - 



Rb \ p{m) ( Rb 



R 



bo 



Rb 



Then T(,,b a £ F( m ) arL d vanishes for r = (1 — rj)Rb , and we consider the differential, Ta = T&.b — P^\ with 
the same domain. The goal is to prove, 



(A.106) 



e mti T A (r)\\ m{R2) <C 



\b-bn 



for some fixed constant C . To do so, consider the equation for Ta written as 

Ri 



L M _ M r 2 i Ta = _ I i I 



i-*v 



fo 2 ,-^ 



b ° ^2 



(A.107) 




^0 



/) 2 



E>2 »0 



—r 2 P (m ' 1 



^b 

bo 



bn ' 



where T^ is the linerized operator from equation (1.10 1. We will use Fbfia to denote the final right hand 
term of (A.107). Note that in the case m — 0, and thus R^ = |, the final multiples of Ta and P^ 
collapse. AH three right hand terms of (A.107) are bounded in the sa me way as in 18 p609], with only 
minor adaptation^] To conclude the argument from 18 and establish (A.106) there only remains to show 
the following Lemma: 

Lemma A.l. Let /i + < be the lowest eigenvalue of L ? , and (f> + G T 2 the corresponding normalized 
eigenvector. For b > sufficiently small with respect to i], and assuming n > is itself sufficiently small, 



L 



(m) 



1 



r ) w,w) > 8 + \\w\\ H i - — (w,(j) + ) , 



for 5+ > constant and any w G H} m \ vanishing at r — (1 — rfjR},. 



^The terms due to i?j> ^ | have no effect. Part of the error term G\(R) that appears in |18| has been moved to the left 
hand side of ( A.107 1, so that the constant Aq that appears in [l8| can be ignored. 
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Lemma A.l is analogous to 19 equation (212)], and is adapted from Lemma 2.15 by using a cutoff and 
the exponential decay of <j>+. Details can be found, [l9j p660]. 
Step 6: Frechet Derivative on Fixed Domain 

The aim is to prove that there exists, 



(A.108) 



W b TbM 



b=bn 



e H (m)- 



We will follow the argument of |T8l p610], and revisit equation (A. 107). In the limit b — > bo we have, with 
respect to L 2 -norm, 



(A.109) 



4 



Note that by direct calculation, 

d 



di 



b-b ( 



b=b 



= 0-3 R 



b-bn 



d 

0, 



b.b t , 



b=bn 




and clearly exists. To show equation (A.108), we recall from Step 5 that, for bo > sufficiently small, 

(m) bl 



L 



is invertible over the subspace of 



functions that vanish at r = (1 — T))Rb- 



Step 7: Uniform Bound for d b T b ^ ba \ b=bg (in C with exponential weight) 

Revisit equation (A. 107 1, again in the limit b — > bo with respect to L 2 norm, 



L 



(m) 



f 2 + 3((^) 2 -K l) ) 2 )) 



db T ' 



b,b„ 



b=bn 



b = bn 



Similar to Step 4, we apply a maximum principle argument on the region - < r < (1 — rj)Rt, to prove, 



||e d-^e(^) ^ 



The full argument is the same as |18 
Step 8: Uniform Bound for dbP, 



(m) 



Let P, 



(A.110) 



(m) 



p( m ^ where 



with the good behaviour, V0f,||£ 



I! < 1 

b=b °o 

p610-611] with only minor adaptations, 
(in C 2 with exponential weight) 

=£>o 

are the smooth cutoff functions, 

1 for r < (1 - rj) 2 R b 
for r > (1 - r])R b , 
— >• 0, as b —> 0. Alternately, 




(A.lll) 

The goal is to prove that, 
(A.112) 



P, 



(m) 



^b - <i>bo)P b 



(m) 



0bo ( 



P. 



(m) 



P. 



(m) 



P, 



('") 



d b P » llc2 



o 



o. 



which is the second precursor to (2.22). Regarding the first right hand term of (A.lll), we may re-express 
in terms of T b bo . Then by Step 7 and the support of cj) b — (f> bo , the contribution from that term is 



(m) 



Pi 

neglible. The remaining term, cf> bo (^P b "^ — P b ^\, is treated with calculations similar to those applied to 

Ta in Steps 5, 6 and 7. Details can be found, [18j p611-612]. 
Step 9: Supercritical Mass 

The proof of (2.23) is d ue to 22 Lemma 1]. Here, we give a summary for the reader's convenience. To 
begin, note from equation (2.18) that is formally a function of b 2 . Then from Step 8 and the chain rule 
we conclude that, with an exponential weight, d^ b ^P^ m ^ is bounded in C 2 . From equation (2.18) it can be 
shown in the limit b — > that, 



(A.113) 



r 9 p(m) _ r p(m) 

+ d{b 2 ) 6 ~ 4 b ■ 
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Consider then a product of | |A.113| by Ai?( m \ 

l - J \x\ 2 \R^\ 2 dx =- l -(r 2 rt m \KR^) 

= -lini(L + a (b2) P b (m) ,A J R(" 1 ) 

This concludes our summary of the proof of Proposition |2.1 

roof 
reads, 



2 / \ 

Proof of Proposition^^ Apply the point transformation, e 



-imO 



r m Z(r). Then equation (2.52) 



where r m * b = A0 h P b (m) + • VP b M 



. The arguments of 19 Appendix E] and 



22 



Appendix A], then prove a version of Proposition 2.3 for e* 6 ^ Z(r), as a radial function on 
accounting for the equivalences of norms, this proves Proposition |2.3| 



p2m+2 



By 
□ 



Appendix B. Details of Numerical Methods 



Our numerical methods closely follow those detailed in 15 , employing the Fortran 90/95 boundary value 
problem software described in 30 . We briefly review it here. 

The software is designed to solve two point boundary value problems of the form 

(B.114) i_ y= Is y + f ( r; v) ; 

dr r 

by nonlinear collocation. Note that the algorithm handles r _1 singuralities. All of our computations were 
performed on the domain [0,50] with tolerance 10~ 10 . 



Codes that can be used to reproduce the computations presented here are available at http://www.math. 
toronto . edu/simpson/f iles/vortex_dist . tgz 



B.l. Point Transformations. Unfortunately, the equation for the vortex state, (1.7), and the operators 
Z'lTad and -C^Tad' hiclude r~ 2 singularities. We address this with the point transformation 



(B.115) 



R(™)(r)=r m M m \r). 



Similarly, U = e imf, r m LK m )(r) for any of the dependent variables. With this transformation, the vortex 
equation becomes, 



(B.l 16a) 
(B.116b) 



(P( 



rn)\// 



2m + l (R (m) )' - R {m) + r 2m (p( m )) 3 



(R (m) Y(0) = 0, lim R( m \r) = 0, 



and the operators £i™ad'^2rad become, 



(B.117a) 



(B.117b) 



£ (m) U - r r ' 



-u" 



= r m dU 



r (m) y _ r, 
i "2.rad Z - r 



-u" 



= r m C 2 Z 



— f/' + 3r 2m R^(mR^ + r(R^)')U 



2n l±lu' + r 2m R^ m \mR^ + r(R^)')U 
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The right hand sides of (3.901 conveniently become, 

(B.118) R (m) =r m^(m) 

(B.119) AR^ = r m {(to + l)R^ + r^" 1 ))'} 

A 2 i? (m) = r rn f ^ + 1 j2^(m) + (3 + 2m )r( J R^) / + r 2 (R^)"\ 

(B.120) J 

= r m { [(to + l) 2 + r 2 ] + 2r(i?( m ))' - r 2(m+l)(£(™))3 I 

B.2. Artificial Boundary Conditions. As the algorithm is designed to compute on finite intervals of 
[a, b], we must compute on [0,r max ], where r max is sufficiently large. This neccessitates the introduction 
of an artificial boundary condition on M m \ the vortex state, and Uj and Zj solving the boundary value 
problems (3.90). The analogous question in the index function computations is verifying that there are no 
zeros beyond r max which might have been missed. 

To develop the artifiical boundary conditions, we examine the asymptotic behaviour of the solutions, using 
that potential terms are exponentially decaying. For the vortex state, 

(B.121) R^ir) ocr~ m ~h- r 

This gives us the boundary condition at r max 

(B.122) 0R (m) )'(r max ) + (l + 2 J^1±1\ = 

which is accurate to 0(r~ 2 x ). 

By similar analysis the solutions to the linear boundary value problems, generically denoted by W, are 

(B.123) W(r) (x r- 2m 

as r — > oo. Thus 

2to - 

(B.124) W'(r max ) + W(r max )=0 

^max 

This too is accurate to 0(r~ ax ). 

B.3. Verification of Results. With these approximations, we solve the following sets of equations, as 
single first order systems: 

• The vortex , and the index functions U and Z, 

• The vortex w m \ the boundary value problem solutions U\ and U2, and the Kj inner products. 

• The vortex RS m \ the boundary value problem solutions Z\ and Z2, and the Jj inner products. 

In computing the index functions, or alternatively the inner products, we are actually solving mixed initial 
value/boundary value problems. 

We now present several a postiori checks on the accuracy of our results. All are based on checking that 
the behaviour of the solutions for large r is consistent with the anticipated asymptotic behavior. 

B.3.1. Verification of the Vortex States. Two related ways of checking that we have adequately computed 
the vortex states are to examine its decay as r becomes large and to see that (B.121) becomes small as 
r — > 00. For the vortices appearing in Figure [2j we plot these two metrics in Figures [4] and [5j With this 
artificial boundary condition, the exponential decay is well captured. 

B.3.2. Verification of the Index Count. In counting the zeros of the index functions from Figure[3J there is the 
concern that there may be another root located beyond r max . To assess this, we note that the asympotically 
free behavior of U and Z is 



(B.125a) U {m ^ ~ c£ m) + C[ m) 

(B.125b) Z< m > ~ + D[ m) 



r -2m 



r -2ro 
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Figure 4. Behavior of the computed vortices as r becomes large. We recover the exponen- 
tial decay. 
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Figure 5. Asymptotically, the artificial boundary condition on the vortex is well satisfied. 

We can estimate these constants by noting 
(B.126a) 

(B.126b) 

(B.126c) 



U 



(m) 



(B.126d) 



-2m 
Z'r 
2m 



D 



(m) 



D, 



(m) 



These constants are plotted in Figure [6] As they show, we have certainly computed into the "free" equation 
zeros m the or functions appearing in Figure | 



regime. More importantly, since C { ™ r > and D { ™ ] < in all cases, we should not expect any additional 



B.3.3. Verification of the Inner Products. For the inner product computations, we verify that in solving the 



boundary value problems, u\ m \ z\ m ^ adequately satisfy the artificial boundary conditions (|B.124[), and that 
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Figure 6. The asymptotical index function constants for winding numbers m — 1, 2, 3. 

the Ki , Ji values are "constant" . The check on the boundary conditions is given in Figures [7] and [8j As 
these figures show, (B.124) is well approximated. 



In computing the inner products, we define 

(B.127) k[ m) (r)= f U^R^rdr. 

Jo 



„(™) i,( m ) Am) -(m) 



"2 ' 3 i 



ji" 1 ' , J2' b \ and jg™"* are defined analogously. Clearly, 



lim k[ m) (r) = K\ 

r—>oo 



(m) 



(B.I28) 

and analogously for the other inner product values. We approximate, 
(B.129) K^n k{ m \r max ), 

for r max sufficiently large that these converge to their limiting values. As Figures [9] and [To] show, this is 
indeed the case. 
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Figure 9. Convergence of the approximate inner product values to their limiting states. 
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